Current relaxation in nonlinear random media 
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We study the current relaxation of a wave packet in a nonlinear random sample coupled to the 
continuum and show that the survival probability decays as P(t) ~ l/t a . For intermediate times 
t < t* , the exponent a satisfies a scaling law a = /(A = x/'oo) where \ i s the nonlinearity strength 
and loo is the localization length of the corresponding random system with \ = 0. For t ^> t* and 
X > Xcr we find a universal decay with a = 2/3 which is a signature of the nonlinearity-induced 
derealization. Experimental evidence should be observable in coupled nonlinear optical waveguides. 
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A fundamental source of physical information arc time- 
resolved decay measurements in quantum mechanical 
systems which are coupled to a continuum via metal- 
lic leads or absorbing boundaries. While the radioactive 
decay is a prominent paradigm, more recent examples in- 
clude atoms in optically generated lattices and billiards 
Q, 0], the ionization of molecular Rydberg states 
photoluminescence spectroscopy of excitation relaxation 
in semiconductor quantum dots and wires 0, and pulse 
propagation studies with electromagnetic waves [j| . Mo- 
tivated by the experimental achievements, the dynamics 
of open quantum systems has also gained considerable in- 
terest from a theoretical perspective and various analyti- 
cal techniques have been developed to s tudy the problem 
in more detail @, S 0, H lH El E3 • One possi- 
ble approach to the problem is to consider the survival 
probability P(t) of a wave packet which is initially local- 
ized inside an open sample of size TV. The total current 
leaking out of the sample is then related to the survival 
probability by J(t) = 

For ballistic/chaotic systems (jj, and for random sys- 
tems in the metallic regime 0, P(t) is by now well un- 
derstood. Recently, also quantum systems with a mixed 
classical phase space have been studied 0, ^| , where 
it was found that P(t) ~ \/t. The same algebraic de- 
cay was found ^3 for disordered (or dynamical) systems 
with exponential localization [lj, 1 1.4 llfi| . In both cases, 
this law is related to localization and tunnelin g e ffects, 
and applies for intermediate asymptotic times [lfj. For 
larger times, localization effects lead to a log-normal de- 
cay of the survival probability ||. 

The subject of the present paper is the survival proba- 
bility in a new setting, namely a class of random systems 
where the evolution is governed by a Discrete Nonlinear 
Schrddinger Equation (DNLSE) (see |l7| and references 
therein), that is 
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i-^Pn{t) = ^ H nm 1p m (t) - xHn{t)\ 2 tpn(t) 



(1) 



be at site n at time t, H nm is a random tight-binding 
Hamiltonian and x 1S the strength of the nonlinearity. 
The nonlinear term in Eq. Q can arise due to a mean 
field approximation for many-body interactions, e.g. in 
the Gross-Pitaevskii framework for Bose-Einstein Con- 
densates or it can result from the description of a 
quantum mechanical particle which moves in a random 
potential and interacts strongly with vibrations |l9| . The 
DNLSE can be also viewed as a description for the en- 
ergy transfer in proteins |2Cj . Finally, in the context of 
optics the DNLSE is capable of describing wave motion 
in coupled nonlinear optical waveguides [2 lL l22t |23| . In 
the latter case the longitudinal space dimension of the 
waveguide plays the role of the time variable. 

Here, we present for the first time the consequences 
on the decay of the survival probability P(t), i.e. on the 
flux out of the sample when coupling a nonlinear ran- 
dom system like Eq. Q to the continuum via absorbing 
boundaries (or conducting leads). In particular, we find 
that for intermediate (but large) times 
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where C is some constant and /(A) is a universal scaling 
function which encodes the interplay between the non- 
linearity and the disorder. We find that /(A < A*) « 1 
while /(A > A*) w 0.35 ± 0.05 where A* - 0(1). 

Moreover, for nonlinearity strengths \ above the de- 
localization border \cr at which localization phenomena 
are destroyed 0, and for times t > t* ~ N 2 /" / (xloo) 2 
(where v w 2/5 is determined by the asymptotic spread- 
ing of the wavepacket's variance var(i) ~ t v ), we find a 
universal asymptotic power-law decay 



p(t)«c(f) 



2/3 



(3) 



Here, \ifj n (t)\ 2 denotes the probability for a particle to 



The intermediate power law (2) is observed for a large 
class of localized initial excitations (i.e. <5-likc, nar- 
row Gaussians) near the surface of the sample which 
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FIG. 1: The survival probability P(t) [Eq. @] shows a power- 
law decay P(t) ~ t~ a for intermediate (but large) times t < 
t*. Shown are three representative examples for k = 5, K = 7, 
N = 1024 and x = 1, 15, 200 (bottom to top). 



have attracted both experimental and theoretical inter- 
est |jl Hoi flit IrH fl3| in recent years. On the other hand 
(3) is independent of the initial condition and applies 
as soon as the wavepacket is spread crgodically over the 
sample. Our results (ll'l.'ili arc confirmed numerically and 
are supported by theoretical arguments. 

To investigate the decay of the survival probability we 
use a modified version [23j of the Kicked Rotor (KR) 
with absor bing boundary conditions [loj |. Based on the 
similarities |l5l between dynamical and Anderson lo- 
calization, it is expected |23j,|2_4j that the same dynamics 
will be generated by Eq. |JT]|, The discrete quantum me- 
chanical evolution from time t to time t + 1 (measured in 
units of a kick period T) is described by the map 

^n(t+l) = ^e^( m - n )j„_ m (fc)e- 4T « m+ « 2 -^^l 2 V m (t) 
m 

(4) 

where the Bessel function J„_ m appears as the re- 
sult of the kick described by the operator [/kick = 
cxp(— ikcos(8)). Here, 9 denotes the angle, n is the an- 
gular momentum while K = kT is the classical kick- 
ing strength. The parameter (f> can be interpreted as an 
Aharonov-Bohm flux through the ring parameterized by 
the coordinate 9. In contrast to the standard KR now 
the kinetic term depends on the wavefunction probability. 
The modified KR J3J models propagation of nonlinear 
waves in optical fibers with a change of the optical density 
inside the waveguide j^,)^. The same model approx- 
imately describes propagation in waveguides with longi- 
tudinal sinusoidal modulation of the boundary [2^, |23| . 

The dynamics of the model for x = is by now 
well studied. The classical motion is chaotic and for suf- 
ficiently large values of K there is diffusion in momentum 
space with diffusion coefficient D ~ ^- ^5|. The most 
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FIG. 2: Power-law exponents a of the survival probability as 
a function of the scaled nonlinearity strength \/\ for various 
localization lengths loo ~ k 2 with a heuristic scaling function 
(full line). Error bars are smaller than the symbol size. Lower 
inset: The scaling parameter A increases linearly with loo (full 
line, slope s=s 1). Upper inset: Power-law exponents a vs. the 
nonlinearity strength x f° r k = 5(+) , 9(A) , 13(D). 



striking consequence of quantization is the suppression 
of this classical diffusion due to quantum-dynamical lo- 
calization |l5|. the dynamical version of the well-known 
Anderson localization 01 . The eigenstates of the as- 
sociated unitary operator U are exponentially localized 
with a localization length ~ k 2 . The dynamical lo- 
calization, however, was found to be destroyed for non- 
linearity strengths \ > Xcr ~ 1/T [23[. A subdiffusive 
growth var(i) « (TxZoo) 475 ^ of the wave packet's vari- 
ance emerged, leading to a uniform spreading over the 
whole sample. The value v = 2/5 was deduced in j2^| 
on the basis of the Chirikov criterion of overlapping res- 
onances and was found to be in reasonable agreement 
with numerical results. For times t > t* ~ p the 

subdiffusion leads to a uniform distribution of an initially 
localized wavepacket over the sample of size N. 

To open the system described by Eq. |j2J we addition- 
ally apply a projection operator V to the wavepacket in 
a fixed interval of momentum states < n < N, i.e. 
\i>(t)) = V\tp(t)). Thus, V describes the complete dele- 
tion (absorption) of the part of the wavepacket which 
propagates outside the given interval. The survival prob- 
ability inside the sample at a time t is then given by 
p(t) = \ip(t)\ 2 = (V>(f)|^>(£)). To suppress fluctuations, 
we concentrate on the geometric average of pit), i.e. 

Pit) = exp((ln(p(i))> ), (5) 

where indicates an averaging over different phases 
(f> (typically more than 20). The initial excitation (un- 
less stated otherwise) is a <5-likc wave packet launched 
at one of the boundaries, i.e. ip n (t = 0) = <5 ni o- In our 
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FIG. 3: Density plots of a wave-packet evolution for (a) 
X = 1 and (b) \ = 50 (parameters as in Fig. 0. Notice 
the development of probability depletion near the absorbing 
boundary for x = 50 which leads to the creation of a potential 
barrier V x = — x\4>n{~t)\ 2 ■ The time axis and the color-coded 
height of the wave function (see inset in a) are both on a 
log-scale. 



numerical calculations, we have used K = kT = 7 and 
k = 3,5, 7, 9, 13 while the sample length N was chosen to 
fulfill always N ^> l^. Due to localization, the survival 
probability shows a decay P{t) ~ 1/t for \ — [53 • 

In Fig. ^ we report the survival probability P(t) for 
k = 5, N = 1024 and three representative values of the 
nonlinear coupling (x = 1, 15, 200). In all cases, P(t) 
clearly displays a power-law decay P(t) ~ t~ a for times 
t < t* . The exponents a for various nonlinearities x an d 
localization lengths Iqo are summarized in Fig. [2] All 
curves have the same functional form, albeit being shifted 
with respect to each other (Fig. [21 upper inset). The 
curves a(x) do coincide, however, when plotting them 
versus A = x/A where A is a scaling parameter (Fig. [21 
main part). We find that the scaling parameter A de- 
pends linearly on (Fig. [21 lower inset) resulting in a 
one-parameter scaling of the power-law exponent accord- 
ing to Eq. 121. Our data indicate that for A < A* ps 0.5 
localization effects are dominant and the survival proba- 
bility decays as P(t) ~ 1/t. In the opposite limit of large 
nonlinearities the decay of P(t) follows a universal law 
P(t) ~ l/t a with a ps 0.35 ± 0.05. We found that the 
smooth transition between the two limits is characterized 
by the scaling function /(A) = 0.35 + 0.65 cxp(-A 3/2 ). 

In Fig. [3 we report the spatio-temporal evolution of 
the wave packet for two limiting values of x in a density 
plot. For X = 1 the initial excitation is essentially local- 
ized in a region close to the absorbing boundary, covering 
the interval < n < in a uniform manner (Fig. Eli). 
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FIG. 4: Decay of P(t) for N = 256, k = 5, K = 7 and 
X = 2, 10, 50 (from top); the numerical result is in agreement 
with Eq. (dashed line). We used a uniform (ergodic-like) 
excitation as initial condition. 



This is the same behavior as found for x = (data not 
shown). As the nonlinearity increases, subdiffusion takes 
over, leading to a spreading of the dominant fraction of 
the wavefunction into the bulk of the sample. This re- 
sults in a gradual depletion of the wavefunction from the 
boundary zone and a creation of an effective potential 
barrier, the height of which increases with the nonlin- 
earity strength x a s V x (t) = — xl^nit)] 2 ■ This barrier 
traps the probability inside the bulk of the sample and 
obstructs the outwards flux thus leading to the observed 
slow decay of P(t) for intermediate times t < t* j^. 

We want to understand the one-parameter scaling of 
the power-law exponent a(xJoo) = /(A) (sec Fig. EJ- 
We recall that the slow decay to the continuum in the 
presence of the nonlinearity is due to tunneling through 
a barrier of height V x ~ xlV'l 2 which is created in a regime 
n ~ loo from the boundary. Assuming that for very small 
nonlinearities the components ipn of the wave packet are 
similar to the case x = 0, i.e. \ip n \ ~ exp(— n/Z 00 )/v / Ic», 
we obtain V x ~ x/^oo- The barrier is totally transparent 
if V x < E where E ~ 0(1) is the energy of the excitation. 
Consequently we get that A* = x* /loo ~ 0(1) indicating 
that the ratio A = x/^oo controls the behavior of P(t) 
and therefore the value of a (sec Eq. ©). 

Next, we give a qualitative argument for the power-law 
decay |j2Jl when A < A* . We start with the perturbative 
expression for the survival probability for x = 01; i- e - 

poo 

P{t) ~ / \c k \ 2 e- T ^dk , (6) 
Jo 

where \c k \ 2 = \(tp(t = 0)\u k )\ 2 - exp(-2fc/Z 00 )/? 00 are 
the overlaps of the initial state \ip(t = 0)) with the ex- 
ponentially localized cigenstates \u k ) of the operator U 
[Eq. (D] for x = 0. The decay P(t) - 1/t results 
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from Eq. © when using 1^ ~ Tk ~ exp(— 2k/l oa )/l oc , 
where Tk is the transmission probability of the k— th 
mode. In the presence of a nonlinearity, it is known 
that the transmission probability still decays exponen- 
tially [24I Er| , and therefore it is reasonable to assume 
that 71. ~ exp(— 2k/l x ) where l x = aloo is an effec- 
tive length determined by both, the tunneling through 
the barrier V x and the disordered potential. Since the 
existence of the barrier additionally hampers the trans- 
mission Tfc with respect to x = 0> it is reasonable to 
assume that a < 1. Substituting this expression for Tk 
into Eq. JBJ we obtain Eq. 

For times t > t* and x > Xcr, the subdiffusion leads to 
an ergodic-like distribution of the initial excitation over 
the whole sample jl^- A simple approach based on a 
modified diffusion equation suggests that P(t) should 
consequently decay in a stretched-cxponcntial manner 
|27| . It is thus surprising that our data contradict this 
very expectation (see Fig.^J. In fact, we find a universal 
power-law decay with exponent a ~ 2/3. This can be un- 
derstood by working in the basis of localized eigenstates 
\uk) of the linear case (x = 0). In this representation, 
we can estimate the transition rate from a quasi-energy 
level Ek to other levels in a discance yj var(t) to be Tk ~ 
(Tx) 2 /[var(0] 3/2 E3 Thinking semi-classically, the lat- 
ter can be associated with the escape time tk ~ 1 of a 
particle that was initially located at a distance -^/var(t) 
from the absorbing boundary. At any time t the number 
[i of particles with inverse escape times Tk > t^ 1 is thus 
given by fj, — ^/var(t) — [(Tx) 2 ^ 3 - As the probabil- 
ity to survive inside the sample up to a time t is given 
by P{t) cx dfi/dt [h|, this leads to P{t) ~ (T X /t) 2/3 in 
agreement with our data (see Fig. . 

We have checked also that for x < Xcr localization 
effects are dominant leading to a decay of P(t) as in the 
case of x = 0,E3 Finally, in the limit £ — > 00, | 2 
becomes very small due to the loss of norm from the 
absorbing boundaries and P(t) decay as for x = [9Hl(ij|. 
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